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Abstract: 

In a first part, we justify the feasibility of substituting a photon leg by 
a neutrino current in the Euler-Heisenberg Lagrangian to obtain an ef- 
fective Lagrangian for the process 7^ — > 77^ and its crossed reactions. 
We establish the link between these processes and the four-photon scat- 
tering in both the Standard Model and the effective theory. As an ap- 
plication, we compute in this effective theory the processes 7^ — ► 772/ 
and 77 — > ^jvv and show how to use the 77 — > 77 results as a check. 
We settle the question of the disagreement between two computations in 
the literature concerning the reaction 77 — * ji/u. In the second part, we 
present results of the direct computation of the photon-neutrino five-leg 
processes in the Standard Model, discuss possible astrophysical implica- 
tions of our results, and provide simple fits to the exact expressions. 
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In a first part, we justify the feasibility of substituting a photon leg by a neutrino current in the Euler-Hcisenbcrg Lagrangian to 
obtain an effective Lagrangian for the process -yu — > 771^ and its crossed reactions. We establish the link between these processes 
and the four-photon scattering in both the Standard Model and the effective theory. As an application, we compute in this effective 
theory the processes yi/ — > 77^ and 77 — > -yvi> and show how to use the 77 — > 77 results as a check. We settle the question of 
the disagreement between two computations in the literature concerning the reaction 77 — > ■yuu. In the second part, we present 
results of the direct computation of the photon-neutrino five-leg processes in the Standard Model, discuss possible astrophysical 
implications of our results, and provide simple fits to the exact expressions. 



1 Introduction 

Processes involving photons and neutrinos are potentially 
of interest in astrophysics and cosmology. However, it 
was realized some time ago inEJ that 4-leg processes (7^ — > 
7^, 77 — > vv and vv — > 77) are too strongly suppressed 
to be of relevance. The reason for this suppression is the 
prohibition of a two-photon coupling to a J = 1 state. 
This is because of Yang's theoremB. On the other hand, 
this theorem does not apply to 5-leg processes involving 
two neutrinos and three photons, such as 
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The extra a in the cross section is compensated by an in- 
terchange of the 1 /Mw suppression by an 1 /m e enhance- 
ment Ba. The relative enhancement of the 5-leg process 
versus the 4-leg one happens to be of several orders of 
magnitude, depending on the energy. 

Iro, Dicus and Repko derived an effective Lagrangian 
for the above five-leg photon-neutrino interactions by us- 
ing the Euler-Heisenberg Lagrangian that describes the 
photon-photon scatteringd. Moreover, in the literature 
there., already existed a computation by Shabalin and 
HieuH of the second pppcess in ([!]) whose result disagrees 
with the one given inu. 

To settle this question, in a recent worku we com- 
puted the first and the second processes (0), in the 
framework of theL effective theory, confirming the results 
reported in ref. In section 2, we justify the feasibil- 
ity of this approacha and give another derivation of the 
5-leg effective vertex starting from the Euler-Heisenberg 
Lagrangiand. 

The effective approach gives reliable results for ener- 
gies below the threshold for e + e~ pair production, while 
its extrapolation to energies above 1 MeV, interesting 
to study, for example, supernova dynamics, is suspect. 



Therefore, an exact calculation of the processes (Q) is 
important to see their role in astrophysics and the range 
of validity of the effective theory. This calculation □ is 
summarized in section 3, with the main results. 

Very recently, parallel work in the same .direction has 
been carried out by Dicus, Kao and Repkota, so that we 
had a chance to compare our numerical results, finding 
complete agreement between the two independent calcu- 
lations. 

We end up with a discussion on some of the impli- 
cations of the exact results of these 5-leg processes in 
astrophysics and cosmology, and we give our conclusions. 

2 Effective theory 

The starting point is the leading term of the Euler- 
Heisenberg Lagrangian El, which describes the photon- 
photon scattering of Fig. la: 
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where F^ v is the photon field-strength tensor and a the 
QED coupling constant. 
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Figure. 1: Four-photon interaction (a) in the effective 
Lagrangian and (b) in the SM. 
In what follows we will show the relation between the 
effective Lagrangian of eq. (|J) and the one describing 
processes (0), which is given by 

£eff = ^[5 (iVi^j (Fx p F x ») - lW pv F^F Xp F^ (3) 
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where 
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and F M „ stands for the field strength of the new "gauge 
field" A v = ^7^(1 — 75)^ = 2r t ,, with T v the neutrino 
current 
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Figure 2. SM leading diagrams contributing to five-leg 
photon-neutrino processes. 
The keypoint of the proof is to demonstrate that the 
SM amplitudes of diagrams (a) and (b) of Fig. 2, which 
we called Aijk and -Byfc, respectively, can be rewritten in 
terms of diagram (b) of Fig. 1. 

The total amplitude M(X,p,cr) of the 5-leg process 
reads 



and l/A z = 1/ (fr 4 +p 5 ) z - Mf J ~ -1/M|. 

The same trick can be applied to the couples of B's. 
Let us take, for instance, 
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The second step of the proof is to shrink diagram 
(b) of Fig. 2 to diagram (b) of Fig. 1 with the photon 
leg substituted by the neutrino current. This is done by 
expanding the W propagator inside the loop, using 
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where A w (q) = (q+k) 2 — and k — p 2 +P3 +P5 . Once 
introduced in eq. ([?]), the first term in the r.h.s. of eq. 
(||) gives rise to an Li type of integral plus a contribution 
of next order in On the other hand, the second 

term on the r.h.s. of eq. (JsJ) is of order and can 

therefore be neglected. 

In conclusion, at leading order in the set of 

four diagrams (from a total of 12) is always proportional 
to Li'. 
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where we denote the photon polarization by A, p and a 
(seeQ for the explicit form of A^k and B^k). 

The reason why the terms in eq. (||) are collected in 
pairs is because, when adding the two terms in each pair, 
and using the reversing invariance of the 7-matrix traces, 
changing q to =q and m e to — m e , the 75 piece of the Zee 
vertex cancels □. For instance, the couple 

A 123 + ^321 = e a (P 1 ,X)e^P 2 ,p)e^(P 3 , a)(A a ^ + A$?), 



where 
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cancels completely its axial part. We are using the nota- 
tions of Refs. 



Q±i = Q=B T m e , 

Q±i = Qo ± Pi , 
D±i = Q ±i ■ Q±i , 
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We have shown that in the SM the 5-leg processes (|lj) 
at leading order correspond up to a global factor to the 
4-photon scattering process; both theories should thus be 
described at this order by the same effective theory, after 
substitution of one polarization vector by the neutrino 
current 
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and fixing the overall constant C. This constant C (eq. 
(^)) is fixed unambiguously by considering the following 
ratios of amplitudes in the large-m e limit 
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where P stands for our 5-leg processes. 



2.1 Cross sections 

Using the obtained effectivaLagrangian eq. @) that co- 
incides with the one used inH we have evaluated the cross 
section for the first two processes: 
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Inta we have also computed the unpolarized differential 
cross sections and the polarized cross sections, using as 
a check for each polarized cross section of the second 
process the well-known polarized cross sections of the 
4-photon scattering process. The results inta coincide 
for the differential and total cross sections with the ones 
obtained inu and are, therefore, in disagreement with the 
ones obtained by Shabalin et al. inQ. 

3 Direct computation of the 5-leg processes 

The range of energy in which the processes (|l|) are rel- 
evant is well below the W mass, so that we treated the 
neutrino-electron coupling as a four-Fermi interaction. 
We also assumed massless neutrinos. The total ampli- 
tude reads 

3 

M(A, p, a) = -2(1 + v e ) r; Y, T i(Pi - P2, Ps, A, p, a) (12) 



large use-af the Kahane-Chisholm manipulations over 7 
matricesl!3. Such identities are strictly four-dimensional, 
while we are, at the same time, using dimensional reg- 
ularization. Our solution is splitting, before any trace 
manipulation, the rt-dimensional integration momentum 
appearing in the traces as □ q — > q + q, where q and q 
are the four-dimensional and e-dimensional components 
(e = n — 4), respectively, so that q ■ q — 0. The 7 algebra 
can then be safely performed in four dimensions, at the 
price of having additional terms. In fact, the splitting 
q — > q + q is equivalent to redefining m 2 — > m 2 — q 2 from 
the beginning. The net effects are then the extra inte- 
grals containing powers of q 2 in the numerator, whenever 



is present in the formulae. The compute 



ion of such 



integrals in the limit e — * is straightforward! . Finally, a 
standard Passarino-Veltman decomposition!^ of the sim- 
ple remaining tensorial structures in terms of scalar loop 
functions, concludes our calculation. We implemented 
the outcoming formulae in a Fortran code, performing 
the phase-space integration by Monte Carlo. Numerical 
results are reported in the next section. 



I2 and l£ are obtained from 1^ by the replacements 
P2 «-> P3 and e(p 2 ,p) «-► e(p 3 ,cr) for 1%, and p 2 <-> p\ 
and e(p2,p) *-* e(pi, A) for 1%. The factor is 
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A detailed description of the method used for the compu- 
tation of the total amplitude eq. ( |l2| ) can be found inQ. 
The reduction of M(A, p, a) to scalar one-loop integrals 
is performed with the help of the technique described 
in El. The general philosophy of such a method is using 
the 7 algebra in the traces to reconstruct the denomina- 
tors appearing in the loop integrals, rather_than making 
a more standard tensorial decomposition E2l. The algo- 
rithm can be iterated in such a way that only scalar and 
rank-one functions appear at the end of the reduction, 
at worst together with higher-rank two-point tensors. In 
that way, all the results are expressed only in terms of 
scalar functions with 3 and 4 denominators, rank-1 inte- 
grals with 3 and 4 denominators, rank-2 integrals with 
3 denominators, and rank-3 functions with 3 denomina- 
tors. This already provides an important simplification 
with respect to the standard decomposition, in that the 
computation of tensors such as 
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is completely avoided. A suitable choice of the polariza- 
tion vectors El as explained inQ is the key ingredient in the 
case at hand. To obtain compact expressions, we made a 



3. 1 Results 

Our formulae remain valid also when including all neu- 
trino species. In this case, only the first diagram in Fig. 
2a contributes, at leading order in u>/m e , because the sec- 
ond one is suppressed by powers of w/m M , r . Therefore, 
the inclusion of all neutrinos can be achieved by simply 
replacing (1 + v e ) with (1 + 3v e ) in eq. (|l2|). However, 
we only considered v e in our numerical results. Effective 
and exact computations are compared in Figs. 3, 4 and 
5, for the three processes. Furthermore, Table 1 shows 
the ratio between exact cross section and a for several 
values of uj/m e , then exhibiting the range of validity of 
the effective theory. 

From the above figures and numbers it is clear that, 
for all three processes, the effective theory is valid only 
when, roughly, u)/m e < 2, as expected. At larger lu, the 
exact computation predicts a softer energy dependence 
with respect to the (ui/m e ) 10 behaviour given by the ef- 
fective Lagrangian. 

All these results are in full agreement with those re- 
ported inH. In order that the results be useful, the only a 
option is to fit the curves in Figs. 3, 4 and 5. The results 
of the fits are 
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2.126291og 2 (r) + 0.4067181og 3 (r) - 0.0298521og 4 (r)], 



a A large-m e expansion has been performed. The first term does not 
give numerical predictions that are useful to extend the effective 
theory beyond w = m e . 
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Table 1: Ratio between exact and effective results for the three 
cross sections. The error on the last digit comes from the phase- 
space Monte Carlo integration. 



cr( 7 7 -> JVV) = ^.7.85491 

ir eff (77 — > ^vv) 
0.3435161og 2 (r) - 0.1140581og 3 (r) + 0.01032191og 4 (r)], 
a(vv -» 777) = r _ 6 . 5737 4 
a cS {vv — > 777) 
0.6898081og 2 (r) + 0.150141og 3 (r) - 0.01233851og 4 (r)] 
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: exp[4.42122 
(r) +0.01032 
: exp[5.27548 - 



where the effective cross sections a eS are given in eq. ( fLl| ) 
and ref. and r = ui/m e . All the above fits are valid in 
the energy range 1.7 < r < 100. 
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Figure 3: "/v — > 771/ cross section in fb as a function of 
u>/m e - 
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Figure 5: vv — > 777 cross section in fb as a function of 
w/m e . 

4 Astrophysical and cosmological possible 
implications 

The 5-leg photon-neutrino processes are of interest in 
astrophysics. The processes 1/7 — + ^77 and vv — ► 777 
can affect the mean free path of neutrinos inside the su- 
pernova core, while the process 77— p 5, ^vv is a possible 
cooling mechanism for hot objects £3. 

Basing their results on the assumption 



(7(71/ — > 77^) = an f ^ 

w " y ' VI MeV/ 



CTo = 10 



" 52 cm 2 
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and on thp_data collected from supernova 1987 A, the 
authors ofEj fitted the exponent 7 in eq. ( |l6| ) to be less 
than 8.4, for u> of the order of a few MeV. The physical 
requirement behind this is that neutrinos of a few MeV 
should immediately leave the supernova, so that their 
mean free path is constrained to be larger than 10 11 cm. 
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The effective theory predicts 7 = 10, while, using 
the curves from the exact calculation, a softer energy 
dependence is observed in the region of interest (see Fig. 
3). A fit to the exact curve gives 7^3 for 1 MeV < u> < 
10 MeV, thus confirming the expectations of 113. 

A second interesting quantity is the range of param- 
eters for which the neutrino mean free path for such re- 
actions is inside the supernova core (of 10 km typiral 
size), therefore affecting its dynamics. Always in it 
was found, with the help of Monte Carlo simulations, 
that for several choices of temperature and chemical po- 
tential, and assuming the validity of the effective theory 
(7 = 10), this happens when ui > 5 MeV. Since the exact 
results are now available, it would be of extreme interest 
to see how the above prediction is affected. More in gen- 
eral, we think that the reactions in (|l|) should be included 
in supernova codes. 

Onlhe basis of the effective theory results, the au- 
thors om, suggested that these processes could also have 
some relevance in cosmology. Consider, in fact, the mean 
number N of neutrino collisions, via the first of processes 
(0) , in_p,n expansion time t equal to the age of the Uni- 
verse li-3: 



N = o~(p/v — > 77J/)n„ c t ,n v — neutrino number density. 

By writing n v and t in terms of the photon energy at 
thermal equilibrium {uj ~ kT), expressed in units of 
10 10 K and denoted by T10, we get 



n u = 1.6 x lO 31 ^ 3 
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* = 277, 



(17) 



When N is large, the neutrinos are in thermal contact 
with matter and radiation, while, for N ~ 1 (namely 
a ~ 10~ 42 T^ 1 cm 2 ), the neutrinos decouple. Applying 



the effective formula in eq. (|llj) , the resulting decoupling 
temperature is Tio ~ 9.5, namely uj ~ 8.2 MeV, there- 
fore outside the range of validity of the effective theory. 
By repeating the same analysis with the exact result, 
we found instead that N becomes of the order of 1 at 
uj ~ 1 GeV, and at these energies, other processes enter 
the game. In conclusion, the five-leg reactions in eq. ([!]) 
are unlikely to be important for a study of the neutrino 
decoupling temperature, contrary to what the effective 
theory seemed to suggest. 



1. D. A. Dicus and W. W. Repko, Phys. Rev. D48 
(1993) 5106. 

2. C. N. Yang, Phys. Rev. 77 (1950) 242; 

M. Gell-Mann, Phys. Rev. Lett. 6 (1961) 70. 

3. D. A. Dicus and W. W. Repko, Phys. Rev. Lett. 
79 (1997) 569. 



4. A. Abada, J. Matias and R. Pittau, jhcp- 
ph/980638^ . 



5. H. Euler, Ann. Phys. 26 (1936) 398; W. Heisen- 
berg and H. Euler, Z. Phys. 98 (1936)714. 

6. N. Van Hieu and E. P. Shabalin, Sov. Phys. JETP 
17 (1963) 681. 



7. A. Abada, J. Matias and R. Pittau, hep 



ph/9808294 



8. D. A. Dicu s, C. Kao and W. W. Repko, |hcp- 
ph/980649S| . 



9. R. Pittau, Comput. Phys. Commun. 104 (1997) 
23 and 111 (1998) 48. 

10. G. Passarino and M. Veltman, Nucl. Phys. B160 
(1979) 151. 

11. F. A. Berends et al., Phys. Lett. B103 (1981) 124; 
P. De Causmaecker at al., Nucl. Phys. B206 

(1982) 53; 

R. Kleiss and W. J. Stirling, Nucl. Phys. B262 
(1985) 235; 

J. F. Gunion and Z. Kunszt, Phys. Lett. B161 
(1985) 333; 

Z. Xu, D. -H. Zhang and L. Chang, Nucl. Phys. 
B291 (1987) 392. 

12. J. Kahane, J. Math. Phys. 9 (1968) 1732; 

J. S. R. Chisholm, Comput. Phys. Commun. 4 
(1972) 205; 

E. R. Caianiello and S. Fubini, Nuovo Cimento 9 
(1952) 1218. 



13. M. Harris, J. Wang and V. L. Teplitz, |astro- 
ph/9707l"l5| . 



14. P. J. E. Peebles, Principles of Physical Cosmology, 
Princeton University Press (1993). 



Acknowledgements 

We thank the authors ofi for having informed us about 
their recent computation. 

We also wish to thank M.B. Gavela, G.F. Giudice and 
O. Pene for helpful remarks. 



References 



5 



